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Abstract 

The U (1) x SU (2) Chern-Simons gauge theory is applied to study the 2-D t—J 
model describing the normal state of underdoped cuprate superconductors. 
The U{\) field produces a flux phase for holons converting them into Dirac- 
like fermions, while the SU(2) field, due to the coupling to holons gives rise to 
a gap for spinons. An effective low-energy action involving holons, spinons and 
a self- generated U{\) gauge field is derived. The Fermi surface and electron 
spectral function obtained are consistent with photoemission experiments. 
The theory predicts a minimal gap proportional to doping concentration. It 
also explains anomalous transport properties including linear T dependence 
of the in-plane resistivity. 

PACS Numbers: 71.27.+a, 71.10.Pm, 71.10.Hf 
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The proximity of superconductivity (SC) to antiferromagnetism (AF) in reference com- 
pounds is a distinct feature of the high-T c superconductors. Upon doping the AF goes away, 
giving rise to SC. At the same time, the Fermi surface (FS) is believed to develop from small 
pockets around (±~ , ±|) fl|], anticipated for a doped Mott insulator, to a large one around 
(71", 7r), expected from the electronic structure calculations. To understand this crossover is 
one of the key issues in resolving the high T c puzzle. For this reason, the underdoped sam- 
ples present particular interest due to the strong interplay of SC with AF. A "spin gap" or 
"pseudogap" has been invoked to explain the reduction of magnetic susceptibility x below 
certain characteristic temperature T* and suppression of the specific heat compared with the 
linear T behavior 0. This gap also shows up in transport properties, neutron scattering, 
and NMR relaxation rate. The angle-resolved photoemission spectroscopy (ARPES) data 
show clear Fermi level crossing in the (0, 0) to (jt, it) direction, but no such crossing was 
detected in the (0,7r) to (tt,tt) direction |§. The observed pseudogap above T c is consis- 
tent with d-wave symmetry. Theoretically there have been two competing approaches: One 
starting from the Mott-Hubbard insulator, advocated by Anderson Q using the concept of 
spin liquid, while the other starting from the Fermi liquid (FL) point of view. Along the first 
line, P.A. Lee and his collaborators have consecutively developed the U(l) || and SU(2) || 
gauge theories, whereas the second approach has been elaborated by Kampt and Schrieffer 
0, and Chubukov and his collaborators ||. 

In this paper we apply to the t — J model the U (1) x £77(2) Chern-Simons bosonization 
scheme for two-dimensional (2D) fermion systems ||. This scheme provides a decomposition 
of the electron field into a product of two "semionic" fields, advocated by R.B. Laughlin ||T0| . 



one carrying the charge (holon) and the other carrying the spin (spinon) . It has been shown 
I [|] that a mean-field treatment of a dimensional reduction of such bosonization procedure, 



m 



keeping the "semionic" nature of spinons and holons, reproduces the exact results obtained 
by Bethe-Ansatz and Luttinger liquid techniques, when applied to the ID t — J model at 
t ^> J. For the underdoped 2D t — J model we neglect the feedback of holon fluctuations 
on the £7(1) field B and spinon fluctuations on the SU(2) field V. The holon field is then 
a fermion and the spinon field a hard-core boson. Within this approximation we show that 
the B field produces a flux-phase for the holons, converting them into Dirac-like fermions; 
the V field, taking into account the feedback of holons produces a gap for the spinons 
vanishing in the zero doping limit, at (±|,±|). A low-energy effective action in terms 
of spinons, holons and a self generated U(l) gauge field is derived. Neglecting the gauge 
fluctuations, the holons are described by a FL with FS given by 4 "half-pockets" centered 
at (±|,±|) and one reproduces the results for the electron spectral function obtained in 
mean field approximation (MFA) flT3| , in qualitative agreement with the ARPES data B 
for underdoped cuprates. Due to coupling to massless holons, gauge fluctuations are not 
confining, but nevertheless yield an attractive interaction between spinons and holons leading 
to a bound state in 2D with electron quantum numbers. The renormalisation effects due 
to gauge fluctuations induce non-FL behaviour for the composite electron, including the 
linear in T resistivity discussed earlier [|TJ|]. This formalism describes a smooth crossover 
upon doping from the long range ordered (LRO) AF state to short ranged (disordered) AF 
state with a gap in the excitation spectrum. The minimal gap is proportional to the doping 
concentration and the gap does not vanish in any direction, 
much smaller 
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The euclidean action for the t — J model in 2D can be represented in terms of a fermionic 
spinless holon field H, coupled to a £7(1) gauge field B, and a spin | complex hard-core 
boson E a , a = 1, 2, satisfying the constraint S*E a = 1, coupled to an SU(2) gauge field V, 
and it is given by |TT|| : 



S = IS drJ2[H* (d T - iB (j) - 6) Hj + zB (j) + (1 - (d T + iV {j)) ap E 
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(ij) 



+ £ [(-tHfe'I™ B H j EUPe J ^ V ) a ^p + h.c)) 

-2S C . S .(B) + S C . S .{V) (1) 



+i (1 - HtHi){l - H*H 3 )(\^UPe S ^ V U^p\ 2 - \ 

where P(-) denotes the path-ordering, 5 the chemical potential for the dopants and 
S C . S .(B) = ^fcPxe^BWBP, S C . S .(V) = ^Tr J d 3 xe^ p (V»d»VP + \V>*V V VP) are the 
Chern-Simons action for the gauge fields (we refer to [[I Ijj for further details). The electron 



i f b , if 



field at site j is decomposed as PJTl[1 : Cj a = e ^ H*(Pe ^ ) a pEjp, where 7 3 - is a straight 
line joining site j to oo in a fixed time plane. 

The (local) gauge invariances of ([!]) are: U(l) : Hj — > Hje lAj , B^(x) — > B^(x) + 
d^A{x), A(x) G R; 31/(2) : -> i^JSj, V^(x) -> i?t(x)y M (a:)i2(a;) + G 
SU(2); and holon/spinon (h/s) gauge: f/j — > Hje^ , — >■ e~ l ^T,j, ^ G R. We gauge-fix 
the //(l) symmetry by imposing Coulomb condition for £?. To retain the bipartite structure 
induced by AF interactions, we gauge-fix the SU(2) symmetry by a "Neel gauge" condition: 

— a x K g)> b'l = + J2- Now we split the integration over the V field into an integration 

over a field V c satisfying the Coulomb condition, Y?L=i d^V^ = (from now on \i = 1,2), 
and its gauge transformations in terms of an SU(2)-vahied scalar field g. Integrating over 
Vq and Bq, we obtain 



V; = £(1 - HJH^gjfg^J^d^Tg (x - j) 
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B^ = B^ + 5B^, SB^x) = ^H^H jd^ig (x - j), 

j 

if B 

where e Ja p = — 1 for every plaquette p, and a a are the Pauli matrices. 



Following the strategy in ID |Tl||, we write down the partition function of holons in a g 
background in terms of first quantized Feynman path integral, and find an a priori upper 
bound on it. We then look for a ho lon-dependent configuration g m , V c (g m ) saturating the 
bound, taken as the starting point to add spinon fluctuations. This can be justified in the 
limit i > J, because the effective mass of holes is very heavy [T5|]. For an arbitrary given 

holon configuration the term (a^gjP(e l ^ i ^ V )gja^)n appears for a fixed link (ij) either in 
the "worldlines" of holons or in the Heisenberg term, but never simultaneously; this permits 
a separate optimization of the two cases (see |TI[] ). Using techniques adapted from the 
proof of diamagnetic inequality [16], assuming translational invariance for the minimizing 



configuration g m , neglecting the quartic pure holon term (5 <C 1) in (ffl), and making use 
of results of fl?|| , it follows that for the optimal configuration (see |18[ for further details): 
V<(g m ) = Mg m ) + J2j ^^arg (x - j)% 
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a 

^-Y)ti:ti 3 - -^—d^axg (x-j] 

j 

and on links belonging to the holon worldlines 



%{9 m ) = E(-l)^^^ t arg (x - j)y (2) 



^fe'^'^l^'j^l, (3) 



while on links in the Heisenberg term 



afgfP^^^Wo^) ~0. (4) 



Here V^(g m ) is the slowly varying part of the SU(2) gauge field related to holons, and 
the physical meaning of (f|) and (§) will be explained later (after eq. (0)). We represent 
gj = exp [-iEi^(-l)'V z arg(j - l)]Rjexp [il{-l)\^a y H*Hj], where R is an SU{2) - val- 
ued field, written in CP 1 form: 

(no summation over j), describing the spinon fluctuations around g™ (for which Rj = I). 
With suitable field redefinition, using the £77(2) invariance of the g measure, the action of 
the t — J model can be exactly rewritten as S = Sh + S s , where 

Sh = Si drJ2H*(d T - (aWRfaR^n - 8)H j 

j 

+ J2HH?e- 1 ^ B+5B H 3 {a^R\(Pe l h,) V+SV )R j af) ll + h.c] 

(ij) 

j 

+ El(l " - H;H 3 ){\(aM(Pe J ^ V+ ™)R*<t$)u\ 2 ~ |}- (6) 

(«> 

Notice that for small hole concentration (5 < 1), V is a slowly varying field. 

We now make the first approximation: suppose in (^j) the fluctuations of the V field, 
due to the spinon fluctuations R are small enough that we can set SV = 0. Since the main 
effect of these fluctuations is to convert the SU(2) gauge invariant spinon field into a semion 
field(see 0), to be consistent, we neglect also the feedback of the holon field on B responsible 
for its semion nature, i.e. set 5B = as well. We believe that the proper account of the 
statistics of gauge-invariant spinon and holon fields is less crucial in 2D than in ID, as we 

will see later. Let us consider the variable R \P(^ V )R 3 = 



ab*^ + a*b* 2 b j2 -ab* a b* 2 + a*b* 2 b* n 
-Oib i2 b ol + a*bnb j2 ab i2 b* j2 + a'&a&ji 



(7) 



where a = exp(| Suj) Vz)- I* 1 t ne hopping term of holons only the diagonal elements of (0) 
appear, a kind of gauge invariant Affleck-Marston (AM) variable [12]; in the Heisenberg 
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term only the off-diagonal elements appear, a kind of gauge invariant resonant valence 
bond (RVB) variable. According to the minimization arguments given above, the mean 
value of the AM gauge variable is s-like, real and close to 1 (see eq.(|)), while the RVB 
order parameter should be rather small (see eq.(§)). using a We now obtain a low-energy 
continuum action for spinons rescaling the model to a lattice spacing e C 1 and neglecting 
higher order terms in e. As it is standard in AF systems we define fij = b* a a a f3bjp and 
T9[: Hj ~ tij + (-ly^eLj, = / < 1, (j)) + More precisely the fields are 



assume L , - _ , , ... , y 

defined on a sublattice, e.g. Qj = j 2 , Lj = L- i+ i - 2l ji = j 2 mod(2) and they describe 

the AF and ferromagnetic fluctuations, respectively. It is useful to write Q in CP 1 form: 
Q = Za(T a pZp, with z a ,a = 1,2, a spin | complex hard-core boson satisfying the constraint 

z a z a /. 



Evaluating the holon contribution in MFA, using the absence of the "0" term in 2D [19 



s.= [ d 3 x-[(d n) 2 + v 2 s {d,nf + n 2 v 2 } - -(^) 2 v 2 , (8) 

J u u 



and integrating out over L, we obtain a low-energy continuum limit NLcr model with action 

±l(d n) 2 + v 2 (d,ti) 2 + n 2 v z 2 }-l 

where g and v s are easily derived in terms of J, t, 5, e. 

To consider the effect of the V field, we replace the NLcr constraint Q 2 = f by a softened 
version adding to the lagrangian a term X(fl 2 — f) 2 , and substitute V 2 by its statistical 
averaging over holon configurations, (V 2 ). Using a sine-Gordon transformation [20|, we 
obtain (V 2 ) ~ —5 In 5 (see |l8j for details). For small J, the coupling constant g is small 



and the system is in the ordered phase; the renormalization group flow in the absence of 
perturbation drives g e ff for large distances towards its critical value; the mass perturbation 
induced by (V 2 ) should then drive the system from the ordered to the disordered phase 
(with short range order only). (One might speculate that if we treat the holons as slowly 
moving impurities, consistent with the known results in the limit t ^> J this would 



lead to a kind of "Anderson localization" considered in ||21|| ). Hence our system should 
exhibit a mass gap m{5) vanishing as 5 — >• 0, and absence of AFLRO (at least for 5 1, 
but sufficiently big). This provides a smooth crossover to the insulating AF regime. The 
subleading perturbation appearing in eq.(U) gives rise to a remnant spin-space uniaxial AF 
interaction in short-ranged AF state. 

We can summarize the above discussion by rewriting the NLa model action in CP 1 form, 
neglecting the short range interactions, as 

S* = f d 3 x- [| (d - A )z a \ 2 + v 2 \(d^ + A^)z a \ 2 + m 2 z* a z a ] (9) 

In the NLcr model without mass term (5 = 0), the constraint z* a z a = f and the symmetry 
breaking condition, e.g. {z\) ^ 0, lead to excitations described by a complex massless field 
S = (z*)z2, with relativistic massless dispersion relation corresponding to the spin waves. In 
the NLa model with mass term the absence of symmetry breaking and the effective softening 
of the constraint lead to excitations described by the spin ^ two-component complex field 
z a , with massive dispersion relations. However, the self-generated gauge field = z^d^z a 
confines the spin | degrees of freedom and the actual excitations are described by a composite 
spin 1 spin-wave field. As we shall see, the coupling to holons will induce deconfinement 
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of spin | excitations. In terms of fields b a , ( slave fermion approach), one realises that the 
z— field in the reduced Brillouin zone with two complex components corresponds (at 5 = 0) 
to appearance of an s + id RVB order parameter with a vanishing gap, in MFA, at four 
points (±f , ±~) ||22|| . The NLa model for spinons can also be derived in this representation. 
The low-energy excitations of this model are fluctuations around these points turned into 
massive ones by the (V^) term. Physically, this is a coexisting ir flux plus s + id RVB 
state. From our estimate we expect the RVB s + id order parameter (^) to be much smaller 
then the AM order parameter (|3]). These features are clearly shown in the numerical MFA 
calculations of 



i f B 

Now turn to holons. We use a U(l) gauge with e being purely imaginary and 

assume spinons are in disordered phase with AM parameter ~ 1. In the rescaled e lattice, 
neglecting higher order terms in e and b, the effective action describes the usual 2-component 



Dirac ("staggered") fermions of the flux phase ||24j| , with vertices of the double-cone disper- 
sion relations in the reduced Brillouin zone centered at (±f, ±f ) (in the e = 1 lattice), with 
chemical potential 5. 

Define the four sublattices: (1) for ji, j 2 even, (2) for j\ odd j% even, (3) for j\ even j 2 
odd, (4) for ji, j 2 odd; they can be grouped into two "Neel sublattices" A = {(1), (4)}, B = 
{(2), (3)}. The holon field restricted to the sublattice (#) is denoted by Set: 

e -*i#(i) + e *i#(4)x 

e -i?#(3) +e if#(2) ) 




7o = a z , 7m = K, a x ),4 = lfl A», $ = ¥*) = 

and assign charge = +l(e^ = —1) to the fields on the A(B) sublattice. Then, neglecting 
short range interactions, we obtain the low-energy continuum action for holons: 

Sh= I d 3 x £ ¥ r \ l0 (d - 5 - e^A ) + t(q - e W^))*W 

J r =l 

Here A^ is nothing but the gauge field for the h/s gauge. We can use S* = S*+Sl to compute, 
as in [P5f| , the gauge field propagator induced by the spinon and holon vacuum polarisation. 
Since the spinon is massive, the corresponding vacuum polarization would be Maxwell-like. 
Hence, in the absence of holons, it would logarithmically confine the spinons. However, 
excitations represented by and describe a FL with a small FS (sp ~ 0(5t)) in the 
reduced Brillouin zone around the points (|,±|). Thus the vacuum polarization exhibits 
the Reizer singularity [EJ and the full gauge interaction is not confining. Nevertheless, since 
we are in 2D, the attractive force mediated by the gauge field is expected to produce bound 
states neutral w.r.t. the h/s gauge, i.e. bound states with quantum numbers of the spin wave 
(for a rigorous discussion of a similar problem, see p6| ) and the electron respectively. 



For this reason, neglecting "semion" nature can be justified to some extent in 2D. Even if 
we neglect the gauge fluctuations, the existence of two bands in the reduced Brillouin zone 
gives rise to a "shadow band" effect (the spectral weight for the part of the pocket facing 
(71", 7r) is greatly reduced) due to the presence of 7— matrices, leading to mixing of fermions 



of these bands. The situation is similar to the slave boson case |T3| 
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To conclude we summarize the main differences between ID and 2D cases. In ID: the 
NLa model of spinons contains a 6 = tt topological term, yielding deconfinement; absence 
of V term makes the spinons massless; there is no B term and there is only one holon band; 
the h/s gauge field A vanishes, hence there is no attractive gauge force between holons 
and spinons, so their statistics appears to be crucial. In contrast, in 2D, without taking 
into account the "semion" nature of holons and spinons, but considering the feedback of 
holons on the SU(2) gauge field, producing the spinon gap, we can already understand 
quite a number of peculiar properties for underdoped cuprates: normal state pseudogap, 
small FS, shadow bands, etc. The main features of the MF calculation [113 survive gauge 
fluctuations. Further consideration of these fluctuations between holons and spinons provides 
a binding force between them, and this composite electron shows non-FL behavior, like linear 
T dependence of resisitivity |TJ[ , and others. More detailed consideration of various physical 
properties within the present model will be given elsewhere 
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